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Strongly correlated quantum systems can exhibit ex- 
otic behavior called topological order which is character- 
ized by non-local correlations that depend on the system 
topology. Such systems can exhibit remarkable phenom- 
ena such as quasi-particles with anyonic statistics and have 
been proposed as candidates for naturally fault-tolerant 
quantum computation. Despite these remarkable proper- 
ties, anyons have never been observed in nature directly. 
Here we describe how to unambiguously detect and char- 
acterize such states in recently proposed spin lattice re- 
alizations using ultra-cold atoms or molecules trapped in 
an optical lattice. We propose an experimentally feasible 
technique to access non-local degrees of freedom by per- 
forming global operations on trapped spins mediated by 
an optical cavity mode. We show how to reliably read and 
write topologically protected quantum memory using an 
atomic or photonic qubit. Furthermore, our technique can 
be used to probe statistics and dynamics of anyonic excita- 
tions. 

By definition, topologically ordered states [ 1 ] cannot be 
distinguished by local observables, i.e. there is no local order 
parameter. They can arise as ground states of certain Hamilto- 
nians which have topological degeneracy and which provide 
robustness against noise and quasi-local perturbations. These 
properties of such systems are attractive for quantum memo- 
ries. However, the local indistinguishability makes measuring 
and manipulating the topological states difficult because they 
are only coupled by global operations. One way to access this 
information is to measure properties of the low lying particle- 
like excitations. In two dimensions, the quasi-particles act 
like punctures in a surface which can have anyonic statistics 
and the topological properties are probed by braiding differ- 
ent particle types around each other. The existence of anyons 
also implies a topological degeneracy 0. Quantum Hall flu- 
ids at certain filling fractions are believed to be topologically 
protected and there is a vigorous experimental effort to verify 
anyonic statistics in these systems Q. A standard approach is 
to perform some kind of interferometry where one looks for 
non-trivial action on the fusion state space upon braiding. This 
is manifested as the evolution of a non-trivial statistical phase 
in the abelian case, or a change in the amplitude of the par- 
ticipating states in the non-abelian case. Some experimental 
evidence consistent with observation of abelian anyonic statis- 
tics in a V = 2/5 filled Quantum Hall state has been reported 
|4| but an unambiguous detection of anyons is still considered 
an open issue 0. 



Spin lattice Hamiltonians can also exhibit topological or- 
der and such Hamiltonians can be built with atoms (6) or 
molecules Q trapped in an optical lattice. A significant 
advantage of using atomic systems is that the microscopic 
physics is well known and there are established techniques 
for coherent control and measurement. Suggestions have been 
made for how one would design anyonic interferometers in 
these systems by using local spin operations to guide excita- 
tion along braiding paths [8 9] \W\ . 

We here present a new approach that directly measures 
topological degeneracy and anyonic statistics using global op- 
erations. The technique involves coupling between a probe 
qubit (single ancilla spin qubit or optical mode) and topo- 
logically ordered atomic spins in an optical lattice. A many 
body interaction between spins is mediated by coupling to a 
common bosonic mode of the radiation field via techniques 
of cavity QED ifTTl fT2l [131 [14 1 or, alternatively, via a com- 
mon phonon mode in ion traps fl5l . Our approach avoids lo- 
calizing and guiding excitations while enabling the measure- 
ment of the statistical phase associated with arbitrary braiding 
paths. 

We also note that recent experiments have demonstrated 
braiding operations on small networks of non-interacting 
qubits encoded in photon polarization fl6l [171 . which gen- 
erates a simulation of anyonic interferometry [18 ) . However, 
since the background Hamiltonian vanishes in such systems, 
they are not protected from noise and the particle interpreta- 
tion of the "excitations" is ambiguous. In contrast, the tech- 
nique developed here allows one to probe directly dynamic 
evolution of anyonic quasi-particles within the parent Hamil- 
tonian. In addition, our mechanism can be used to perform 
reading and writing of qubits initially encoded in light or 
atoms into topological memory, which may be useful for of- 
fline storage during a computation and for applications in long 
distance quantum communication (T9ll20ll . 



ATOMIC AND MOLECULAR SPIN LATTICES IN OPTICAL 
CAVITIES 

We focus on physical systems in which a two-dimensional 
optical lattice is placed within a high-finesse optical cavity 
as illustrated in Fig. [2^. To be specific, we consider the 2D 
square lattice Hamiltonian introduced by Kitaev [21 1 where 
each edge of the lattice represents a spin- 1/2 particle (see 
Fig. [T^). Each vertex v or each face / is associated with an 
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operator H v = n ;Gstar ( v )Oj or Hf — Ylj^f&j. These operators 
collide on an even number of edges and hence mutually com- 
mute. We seek to encode in the +1 coeigenspace of these lo- 
cal stabilizers by choosing the so-called surface-code Hamil- 
tonian: 
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(U,J > 0). The ground states of H sm { have a degeneracy 
dimi^,, = 2 2g+h where g is the genus of the surface and h 
is the number of holes 11221 . Designing lattices with genus 
g > 0, such as the surface of a torus, is challenging, but it is 
possible to create several holes (h > 0) in a planar lattice by, 
for instance, deactivating regions of the lattice with focused 
far detuned lasers. Alternatively, the planar code with spe- 
cific boundary as shown in Fig. [1J provides a ground state 
degeneracy of 2. The logical states are coupled by the opera- 
tors: Z = Ylj e c z a Z j andX = Ylj e( ^ x c X j where the configurations 
Cz(Cx) are strings on the lattice (dual lattice) as illustrated in 
Fig. [1}!. 

There are several experimental proposals to implement the 
spin lattice Hamiltonians with topological order. For ex- 
ample, Kitaev's honeycomb lattice Hamiltonian Z/hcb ( see 
Fig. [Tp)[23| can be designed in optical lattices with ultra- 
cold atoms using controlled spin exchange interactions (6), 
or with molecules using microwave induced dipole-dipole in- 
teractions Q. With an appropriate choice of coupling param- 
eters ll23ll . the honeycomb lattice Hamiltonian has a gapped 
abelian phase with a low energy effective Hamiltonian locally 
equivalent to H sur f. In the following, we will assume the sys- 
tem interacts via H sur f, but our results are also applicable to 
other spin lattice Hamiltonians. 

We now consider how to implement the global operations 
for the spin lattice system. In particular, we are interested in a 
specific type of global operation: products of Pauli operators 
on a set of spins whose corresponding edges in the lattice form 
a connected string. Such global operators are called string op- 
erators. For example, the generators for the encoded qubits 
(Z and X) are string operators (see Fig.[T^). All string oper- 
ators are equivalent to S z c = Yljec a 'j U P t0 local single spin 
rotations, where C is the set of selected spins. For example, 
5^. = TljecGj = Y[jecHjO Z jHj where Hj is the Hadamard ro- 
tation for the jth spin. 

In our setup, the topological memory consists of a spin lat- 
tice of trapped atoms or molecules inside an optical cavity as 
illustrated in Fig. [2^. The off-resonant interaction between 
the common cavity mode and selected spins is described by 
the quantum non-demolition (QND) Hamiltonian E4ll25l : 



H 



l^aY^O), 



(2) 



which is achieved by choosing the cavity mode with a large 
detuning A from a spin-dependent optical transition as shown 
in Fig. ph. The coupling strength is % — g 2 /2A where g is 
the single-photon Rabi frequency for the cavity mode. The 




FIG. 1: Generators for the encoded qubits. (a) A planar code which 
encodes one logical qubit in the ground states. There is a spin- 1/2 
particle (filled dot) for each edge of the lattice. The interactions of 
the local Hamiltonian H sm f are along edges that bound a face /, and 
edges that meet at a vertex v. The strings Cx.z indicate paths of 
products of CT 1 ' 2 operators that are logical operators on the code, (b) 
A nearest neighbor local Hamiltonian Z/hcb on the honeycomb lattice, 
with a spin- 1/2 particle for each lattice site. The (green, red, blue) 
edges represent interactions of type (o^cr 1 ,a y a y ,<7 z O z ). In the limit 
that the interactions along the blue links are much stronger than those 
along the other links, the ground subspace has a gapped Z2 topolog- 
ical phase 1231 . Physical (<T r ,o y ,O z ) spin operations as part of the 
strings Cx.z are indicated by bold (green, red, blue) circles around 
the spins. At qubit k, the string crossing, the operation is CT?ojL 



QND Hamiltonian preserves the photon number n a = a* a of 
the cavity mode. In addition, the cavity mode also interacts 
with an ancilla spin that probes anyonic statistics associated 
with quasi-particles. 

Similar to the previous schemes |[8][l0j to measure anyonic 
statistics, we assume selective addressing of spins in the lat- 
tice so that we can perform single spin rotations as well as 
switch on/off the coupling between the cavity mode and the 
spins. Such selective manipulation can be achieved using ad- 
dressing lasers with shaped intensity profiles GcjI [271 . The 
key new ingredient, however, is that we use the common cav- 
ity mode to mediate global string operators. In this way, we 
avoid problems involving maintaining adiabaticity and local- 
ization while braiding quasi-particles. And most importantly, 
we are able to achieve controlled-string operations ( A [S^ z ] ) 
for an arbitrary string C. 

The idea of controlled-string operations can be illustrated 
by considering a situation when the cavity mode is first pre- 
pared in some superposition of zero and one photon states. 
Within this subspace, the evolution of the QND Hamiltonian 
for interaction time i = ti/2% yields 



U = ex V [-iHx}= Hf c rK' C3) 
I for n u = 

(-0 Wc n>6co5 for« a = l 

where Nq is the number of elements in C, and the second 
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FIG. 2: (a) Inside a cavity, an optical lattice carries spins for topo- 
logical memory, with individual spin addressability |26 27], (b) 
The energy levels of a selected memory spin interacting dispersively 
with the cavity mode, which implements the QND Hamiltonian of 
equation ([2}. The coupling coefficient is % = g 2 /A, with single- 
photon Rabi frequency g and detuning A. (c) The energy levels 
of the ancilla spin and the cavity mode for the single photon ap- 
proach. The control laser fl (f ) connects the states | 1) A ® |vac) and 
I 1')a ® a ^ l vac )> an d enables coherent creation and absorption of a 
cavity photon conditioned on the ancilla spin, (d) Cartoon illustration 
of the procedure for the implementation of single-photon approach 
for controlled-string operations: (1) Initialize the ancilla spin (the 
left highlighted spin) in a superposition state a |0) A + (3 1 1) A (blue for 
|0) A and red for \ l) A ), with no photon in the cavity and state 
for the topological memory. (2) Coherently create a cavity photon 
(orange shade) for ancilla spin state | 1) A (upper branch); no photon 
is created for ancilla spin state |0) A (lower branch). (3) Switch on 
the interaction between the cavity photon and the selected spins. If 
there is a cavity photon (orange shade), a non-trivial evolution 
(pink dots) is implemented. (4) Turn off the interaction and coher- 
ently absorb the cavity photon into the ancilla spin. Finally the state 
OC |0) A ® \\\f) s + p | l) A ® S z c |\|/) s is prepared. 
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equality uses the identity exp 

tary evolution will apply the string operator to the topo- 
logical memory, conditioned on one cavity photon. With 
such controlled-string operations, we can conveniently access 
the topological memory, and build anyonic interferometry to 
probe braiding statistics and dynamics of quasi-particles. 

In practice, however, it is actually easier to control the an- 
cilla spin rather than to directly manipulate the photon number 
state. Therefore, in the following, we will present two ap- 
proaches to controlled-string operations between the ancilla 
spin and the topological memory. 



using spectroscopically resolvable energy levels as shown in 
Fig. |2j;. Starting with no photon in the cavity mode |vac) 
and ancilla spin in state a|0) A + p|l) A , we can coherently 
couple the number state of the cavity mode with the state 
of the ancilla spin by adiabatically increasing the Rabi fre- 
quency £2 (t ) of the control laser until it is much larger than 
the single-photon Rabi frequency g' . The intermediate state 
is then a \0) A <g> |vac) — P | \') A <g) a 1 |vac), having the photon 
number fully correlated with the ancilla spin. Applying the 
QND interaction with the intermediate state realizes the de- 
sired controlled-string operation conditioned on the state of 
the ancilla spin. Finally, we can reverse the state mapping by 
adiabatically decreasing the Rabi frequency, which coherently 
annihilates the photon of the cavity mode and restores the an- 
cilla spin to its logical subspace spanned by {|0) A , | i) A }- Fol- 
lowing the procedure summarized in Fig. [2]J, we can achieve 
the controlled-string operation: 



|1) A (1|®^ + |0) A (0|®I. 



(4) 



The second approach to controlled-string operations is 
based on the idea of geometric phase gates [28 1. Here, the 
bosonic field of the cavity mode starts in a coherent state, 
rather than a superposition of zero and one photon states. 
If our transformation restores the bosonic field to the initial 
coherent state, the entire system will accumulate a quantum 
phase (geometric phase), which is twice the area enclosed by 
the trajectory in phase space of the bosonic field. We activate 
the geometric phase gate using an ancilla spin which experi- 
ences the QND interaction with the cavity mode that can be 
selectively turned on and off l26ll27l . As illustrated in Fig. [5] 
and detailed in the Methods section: if the ancilla spin is in 
state |0) A , the enclosed area vanishes; if the ancilla spin is 
in state |1) A , the enclose area has a different sign depending 
on whether the topological memory is in +1 or —1 subspace 
associated with the string operator S z c , yielding again equa- 
tion 

Various imperfections such as photon loss and deviation of 
the QND interaction can degrade the controlled-string opera- 
tion. However, we can use a cavity with high Purcell factor P 
to suppress photon loss ll29l [30l . and apply quantum control 
techniques to suppress the deviation of the QND interaction 
to arbitrarily high order ||3Tl[32l . As derived in the Methods 
section, the error probability for controlled-string operation is 
approximately yJWcfP, with Nq for the length of the string. 



ACCESSING TOPOLOGICAL QUANTUM MEMORY 



CONTROLLED-STRING OPERATIONS 

The key operation of the single photon approach is the evo- 
lution of the QND interaction described by equation ([3J. In 
addition, the cavity mode interacts with a single ancilla spin 



Controlled-string operations provide an interface between 
the probe qubit which features easy access and efficient ma- 
nipulation, and the topological memory which provides good 
storage. To store quantum states we require two operations: 
the SWAPj n gate which swaps the state of a probe qubit A to 
memory M initialized in |5) M and the SWAP 0Ut which swaps 
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back to a probe qubit prepared in \Q) A . 

SWAP ia =H A -A[z]-H A -A[x], 
SWAP out =A[l] ■H A -A[Z)-H A , 

where H A is the Hadamard gate acting on the probe qubit, 
and A \§\ represents a controlled-string operation. In addi- 
tion, universal rotations of the encoded qubit (generally, ar- 
bitrary unitaries generated by string operators) over the topo- 
logical memory can be implemented either by teleportation 
of quantum gates or by direct geometric phase gate. (See 
SOM for details.) We remark that the ancilla spin can also 
be entangled with another ancilla spin from a different cav- 
ity via probabilistic entanglement generation, and therefore 
our topological memories can be used for quantum networks 
(331 [331 [35] (361 [371 . 

To compare the topological memory and unprotected 
single-spin memory, we introduce the decoherence rate q for 
the unprotected spin due to low-frequency noise. The topo- 
logical memory can significantly reduce the decoherence rate 
by a factor of (8h/J) N <C 1, where S/z is the magnitude of 
the noise perturbation on individual spins and N is the length 
of the minimal string associated with the generators for en- 
coded qubits [21 1. Meanwhile, since the topological memory 
is not protected from errors associated with controlled-string 
operations (with probability ~ y^N/P), we should take them 
into account. Therefore, in terms of total error probability, 
the topological memory outperforms the single-spin memory 
for storage time t > ~ y^N/P. (See the Methods section for 
details). 



ANYONIC INTERFEROMETRY 

We now describe how to use controlled-string operations 
to extract the statistical phase acquired when braiding abelian 
anyons. The definition of anyonic statistics usually relies on 
the adiabatic transport of quasi-particles around each other 
ll38l . with the required condition of adiabaticity to keep the 
system in the same manifold of excited states and prevent ex- 
citing additional degrees of freedom. Note that this procedure 
relies explicitly on the existence of the Hamiltonian. This is 
fundamentally different from anyonic simulation approaches 
|[T6l [T71 [TSl not using topological Hamiltonian, which only 
probe the non-trivial commutation relations of spin opera- 
tors and initially entangled quantum states. However, anyonic 
statistics is a property of quasi-particles associated with the 
Hamiltonian and not just with some specially prepared initial 
state. 

For our spin lattice system with // sur f, there are two types 
of anyons ETl : (1) z-particles that live on the vertices of the 
lattice and (2) x-particles that live on the face (see Fig.|4^,b). 
These anyons are created in pairs (of the same type) by string 
operators: = 5f \%) and |\|f*(/)) = S} \%), where \%) is 

some ground state of the spins, and S| = n ; - e /<7^ and S*, = 
Ylj e it&j are string operators associated with string I on the 



(a) 

4. 



(c) 



MA 















ia 


+ |ocp 


-ia 






< 







(b) 







-ia 


ia 




p 



FIG. 3: Phase accumulation for the approach with geometric phase 
gate [equation |7|]. We use and \l) s to represent +1 and —1 
subspaces of memory spins associated with the string operator S : c , 
respectively. (a)(b) When the ancilla spin is in \0) A state, the en- 
closed area vanishes. When the ancilla spin is in 1 1)^ state, (c) for 
the subspace the enclose area is |cq3|; (d) for the subspace \l) s 
the enclosed area is — 1 0c(3| . The quantum phase accumulated is twice 
the area enclosed. 



lattice and string I' on the dual lattice, respectively (see Fig.|4j>. 
In our approach, string operators can be used to effectively 
move quasi-particles quickly along the string trajectory but 
without exciting other quasi-particles. For example, effective 
motion of quasi-particles with/without braiding is shown in 
Fig.|4^,b. This evolution is described by 



S x l tUt i S Z [JJt 1 S x l iUt\S'i i ^initial) — e ' tc ' 



^initial} 



(5) 



where we introduce time delays, represented by unitary evo- 
lution U t , between string operations. The goal of these de- 
lays is to check that the system stays in the manifold with 
a fixed number of quasi-particles where time delays lead to 
only a trivial dynamical phase. On the other hand, if the 
string operator were to create some complicated intermedi- 
ate states, time delays would lead to complete decoherence. 
The total phase e' Qtot includes both the dynamical contribution 
_ g iZ/(fi+2f 2 +f 3 ) and the stat j st i ca i contribution e' e = -1 

(or +1) in the presence (or absence) of braiding. Therefore, 
we can unambiguously measure the statistical phase if we can 
measure e' 9 "" for both cases. 

The following interference experiment can be used to mea- 
sure the phase e' e ' ot . First, we prepare the probe qubit in a 
superposition state (|0) + |l))/\/2. We then use controlled- 
string operations to achieve interference of the following two 
possible evolutions: if the probe qubit is in state |0), no op- 
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FIG. 4: Braiding operations. (a)/(b) We can move x-particles and 
z-particles in tangled/untangled loops using string operators to im- 
plement operations with/without braiding of anyons. (c) We can also 
apply generators for the encoded qubits to achieve the braiding op- 
eration. The braiding statistics of anyons should be invariant under 
non-crossing deformations of the loops |39|. 




Time t 



eration is applied to the memory spins; if the probe qubit 
is in state |1), the operation S*, U t3 Sj U^S^U^Sj is applied 

to the topological memory, which picks up the extra phase 
giBtot we wan t to measure. After the controlled-string oper- 
ations, the probe qubit will be in state (|0) +e' 6tot |1))/V5. 
Finally, we project the probe qubit to the basis of \£,±) = 
(|0) ±e'* |l))/\/2 with (j) G [0,2%), and measure the opera- 
tor Gty = (^+| — (£- 1. The measurement of (g§) v.s. 
(j) should have fringes with perfect contrast and a maximum 
shifted by (j) = 8 to t. In fact this scheme can be used to mea- 
sure abelian statistics for an arbitrary finite cyclic group as 
described in the Methods section. 

It is crucial to verify that the outcome of the anyonic in- 
terferometry is invariant under repeated experiments with de- 
formed string operators (see Fig. [4]) [39). For example, the 
two ground states of the 2D compass model [40] are coupled 
by perpendicular global X and Z string operators and the phase 
measured using the interferometry scheme above would also 
yield a phase — 1 due to the anti-commutation relations at the 
crossing spin. Yet the ground states are not topologically or- 
dered because deformed string operators do not preserve the 
ground subspace. Since our anyonic interferometry can test 
all possible braiding operations, we can unambiguously ver- 
ify the topological property of anyons. 

Various imperfections will degrade the signature of anyonic 
statistics. The string operators may have errors that excite un- 
wanted anyons, and the topological memory may not fully re- 
store to the ground state after braiding. In addition, the topo- 
logical memory may have anyons from imperfect initializa- 
tion. If these anyons are enclosed by the braiding loops, they 
will affect the phase factor associated with braiding. However, 
neither of these imperfections will prevent us from probing 
anyonic statistics, since they only reduce the contrast of the 
anyonic interferometry without shifting the fringes of (cty). 
We may even distinguish the two types of imperfections from 
the contrast that depends on different loops. The reduction of 
the contrast is proportional to the length of the loops for er- 
rors from string operators, while it is proportional to the area 
enclosed by the loops for errors from imperfect initialization 



FIG. 5: Fringe contrast of anyonic interferometry as a function of the 
time for anyonic diffusion. The fringe contrast quickly reduces due 
to anyonic diffusion (black dotted line). However, we can extend 
the fringe contrast to longer times by applying one (red dotdashed 
line), four (green dashed line), or ten (blue solid line) pairs of time- 
reversal operations of U% within the time interval T. For clarity, we 
only consider the diffusion of two intermediate x-particles induced 
by the perturbation It*®? from equation l|6j. We choose the unit of 

time so that the perturbation strength is normalized J (^(hg) 2 ^ = 1 , 

and we assume that the random field h x e is a Gaussian random process 
with correlation time x c = 10. (See SOM for details.) 



(see discussion in the Methods section). 



PROBING AND CONTROL ANYONIC DYNAMICS 

Our anyonic interferometry provides a tool to study the 
dynamics of anyons. First, consider repeating the pro- 
tocol [equation (|5]l] for anyonic interferometry with the 
time delays {*/} =12 3 between the four controlled-string 
operations. Processes of anyonic creation, propagation, 
braiding, and annihilation will induce a time dependence 
of the final state wavefunction in a general expression: 
\V f im i) = a({tj})\W initial ) + ^({tj})\W ± ({tj})), where 

initial]^ ±{{tj})) = 0- Since the reduced density matrix 

if 1 «({';}) A 
of the probe qubit is p = j ^ a*({f }) j )> we can 

measure the complex coefficient a({tj}) using quantum state 
tomography [41 ] of the probe qubit. 

Probing anyonic statistics can be regarded as special cases, 
with a = e'( e+r l) or e ,T| , and p = 0. Although the anyons are 
immobile for the surface code Hamiltonian, the mobility of 
quasi-particles may change when we include local perturba- 
tions, because the excited states with anyons are highly de- 
generate and any small perturbation to the Hamiltonian can 
dramatically change the eigen-wavefunctions. Consider for 
example a specific diffusion model for anyons induced by the 
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local perturbation 



Hpert — 



E E 

ae{x,y,z} eGAll spins 



(6) 



where the /ig' z field components cause diffusion of x(z)- 
particles and the hi field component causes diffusion of dy- 
onic particles (pairs of neighboring x-particle and z-particle). 
(See SOM for details). The nature of the perturbation (e.g., 
time independent or changing with time) determines diffusion 
dynamics of anyons, which can be observed from the coeffi- 
cient a({tj}) using our anyonic interferometry. 

In addition, we can even control the diffusion dynamics 
of anyons. We introduce the effective time-reversal oper- 
ations U* = ILeAii spins ° z e and U% = UeeMl spins <%> which 
anti-commute with h x e <5 x e and h z e a z e terms of H pert , re- 
spectively. The combination of these operations (e.g., 

x/4 x/4 x/4 t/4 

U£ U£ U\ U£ ) is analogous to spin-echo 

pulses in NMR, which can effectively reverse anyonic dif- 
fusion caused by static perturbations and consequently ex- 
tend the fringe contrast of anyonic interferometry to longer 
time delays, as illustrated in Fig. [5] In essence, by applying 
these global operations, we can localize the anyons without 
any trapping potential. (See SOM for details.) Note that the 
anyonic interferometry is closely related to the Ramsey exper- 
iments in atomic physics |42|, which can now be performed 
with anyonic quasi-particles. 



interesting to consider adapting the present protocol to spin- 
lattice systems with non-abelian anyons |47|. 



METHODS 
Derivation of the geometric phase gate 

We describe the necessary elements to construct the geo- 
metric phase gate illustrated in Fig. [5] First we require the 
displacement operator D(^) = e& ~^* a that can be obtained by 
injecting coherent states through cavity mirrors. The ampli- 
tude and phase of the injected field determine the phase space 
displacement of the bosonic field by 

Second, we need the displacement operation that depends 
on the state of the memory spins: 



„(. 



I) ia if(S z c )=+l 
D(-ia) if (Si.) = -1 



where we use (5^) = ±1 to represent the ±1 subspaces of 
the memory spins associated with the operator S z c . We can 

achieve D (ae'? s ~c^J by applying the QND Hamiltonian for 

time tc — Tt/2% both before and after the displacement opera- 
tion D (ae !1 M. The justification is based on the identity 



D{ae 



R{QO)D(ae ii >)R(-QO), 



OUTLOOK 

Controlled-string operations can be applied to other lat- 
tice Hamiltonians as well 14311441 . which may provide robust 
quantum memory with their degenerate ground states. For ex- 
ample subsystem codes ll43ll can be constructed out of 2D and 
3D nearest neighbor spin- 1/2 interactions that are realizable 
with atomic systems J6] |7|. Our approach can be adapted to 
perform the logical operations generated by strings or planes 
of Pauli operators in the 2D and 3D subsystem codes, respec- 
tively. In addition, the ability to measure global operators on a 
spin lattice provides a means to probe other properties of topo- 
logical phases. For example, a class of topologically ordered 
spin states known as string net states [45 1, which includes the 
ground states of H sm i , have the property that they are invariant 
under large closed loop operations. In the present case, these 
operators are -Xioop(Zioop) = Y\.j^c dmei ° X Z wn i cn have expec- 
tation value 1 . A perturbation on H sur f in the regime U ^s> J, 
by, e.g. a magnetic field, acts like a string tension that reduces 
the amplitude of large loops (on a vacuum reference state). 
In fact there are two phases as a function of the strength of 
the perturbation. For very weak perturbations its has been ar- 
gued that the loop order parameter scales with the perimeter of 
the loop while for strong perturbations it scales with the area 
||451 . These are known as deconfined and confined phases in 
analogy to lattice gauge theory and are examples of phenom- 
ena that could be observed using our protocol. It may also be 



with R(x) 
[0,a] = 0. 

Ji>+iQO 



e ixa'a an( j j-jjg two commu ting operators 
For O = L/ec°; and 8 = 7c/2, we have 



e' 9 n, e cexp 



i-a z 

2°; 



M-rm.IT) n 7GC o}, With m = 
Nc- Therefore, by choosing (j) = (m + l)jt/2, we obtain 

Third, we need dispersive coupling between the bosonic 
field and the ancilla spin (with two levels {\Q) A , 1 1) A }) 

V A = %Aa f a\l) A {l\, 

with coupling strength %a, which can be switched on and off 
via optical control |26, 27 1 or mechanical displacement of the 
ancilla spin. With such dispersive interaction, we are able to 
obtain the displacement operation conditioned on the state of 
the ancilla spin, D (|3 1 l) A (1|) = |0) A (0| <g>I+ 1 1) A (1 1 <8> D (P), 
by the following procedure: (1) apply the interaction V A for 
time t A — k/% a , (2) displace the bosonic field by — p/2, (3) 
apply the interaction V A again for time t A , and (4) displace the 
bosonic field by p/2. The steps (1-3) displace the bosonic 
field by =fP/2 for the ancilla spin in state \Q) A and \l) A , re- 
spectively. Combined with the displacement p/2 from step 
(4), we have the operation D (P 1 1 ) A ( 1 1 ) . 

Finally, the controlled-string operation is a combination of 
the above elements: 

(7) 



7 



The bosonic field is restored to its initial state, while accumu- 
lating a phase depending on both the state of the ancilla spin 
and the value for the string operator as illustrated in Fig. [3] 



Fidelity of controlled-string operations and topological memory 

To evaluate the advantage of using topological memory 
storage, we compare the improvement of storing a qubit in 
spins prepared in the ground states of H sm f (assumed at zero 
temperature) versus the decoherence rate for encoding a qubit 
in a single spin. For long-time storage of quantum memory, 
we expect to gain from the robustness of the topological mem- 
ory can significantly reduce the decoherence rate by a factor of 
(5h/J) N <C 1, where S/i is the magnitude of the noise pertur- 
bation on individual spins and N is the length of the minimal 
string associated with the generators for encoded qubits ETI . 
The controlled-string operations that are used to swap quan- 
tum information in and out have errors including cavity decay, 
radiative decay, and non uniform dispersive shifts in the QND 
interaction. 



and selected spins, described by the following perturbation 
§H = %a f a 52 S/ a /j 

where bj is the relative deviation for the jth spin. In the 
presence of cavity excitation, the implementation of the gate 

Uj = exp I /8cr ) on the jth spin could lead to the gate Uj = 



exp 



i (1+8/) Bo} 



ej = Wj-Uj 



. We define the error by the operator norm 

9 1 8/ 1 ll48l . Thus we have the error for the 
QND interaction, P QND = \\UjecUj -UjecUj\\ < 

This kind of error due to inhomogeneity can be compen- 
sated in two ways. First one could add an optical potential to 
the system which is shaped to equalize the couplings 8 7 for all 
the spins along the configuration path C. A second option is 
to use composite pulse sequences on the system in which case 
it has been shown that the error can be reduced to O (|8y| ) 
for V; using 0(k 3 ) pulses (32) ■ Therefore, P QND ~ N C Q |8|* 
effectively suppressed. 

Summary 



IS 



Errors due to photon loss 

The photon loss is attributed to two physical processes: the 
spontaneous decay with rate y for the optically excited state 
\e) (Fig. [Z^), and the cavity loss with rate K during the QND 
interaction. For single photon approach, the interaction time 

is x = k/2% and the effective spontaneous decay rate is sup- 

2 

pressed to £jy by having large detuning A 3> g for each se- 
lected spin (Fig. |2j)). Under the assumption that the selected 
spins decay independently, the total probability for photon 
loss is 

kx + A?|2 yx > 2k^N/P = P, oss , 

where we define the Purcell factor P = g 2 /KJ f29] [30]|. And 
the minimum probability can be achieved by choosing optimal 
detuning A = gy/Ny/K. 

For geometric phase gate approach, we can choose |a| 2 = 
I p| = Jt/2 so that the total probability for photon loss is 
Pioss ~ |oc| 2 (kx + A^YX^ with x « KA/g 2 . Similar to sin- 
gle photon approach, the probability Pi oss can be significantly 
reduced by having a large Purcell factor P > N. For coherent 
states, we cannot identify photon loss events unambiguously, 
but we can still characterize the errors associated with the pho- 
ton loss. 



The deviation of the QND interaction 

The dominant deviation of the QND interaction is from the 
fluctuations of the coupling strength between the cavity mode 



Combining all the decoherence mechanisms, the error prob- 
ability for the swap in- swap-out process with memory storage 
time t is 

Ptopo-mem = (hh/jf qt + 4ky/WJP + Nt, 

with the pre-factor X = 2k (and 4k 2 ) for the single photon (and 
geometric phase gate) approach. Compared with the storage 
error without topological encoding p K f- mem =qxt, for stor- 



age time t > 



, our topological memory 



9 i 
outperforms storage in a single quantum system. 

We remark that for single photon approach, photon loss in- 
duces leakage errors that can be detected without compromis- 
ing the state stored in the topological memory. Such detected 
errors can be overcome by repetition, which is applicable 
to probabilistic operations such as entanglement generations 
in quantum repeater [19] and distributed quantum computer 
11361 . Detected errors can have a very high tolerable threshold 
for deterministic quantum computation schemes [49|. 



Fringe contrast of the interferometer in the presence of 
excitations 



We refer to anyons left from the initialization as quenched 
anyons, which can result in measurable effects to the phase 
measurement associated with braiding [equation Q]. To be 
specific, we will consider the planar code, and assume that the 
probability to have one pair of initial anyons is p while ne- 
glecting the case with multiple pairs of anyons. If the anyons 
are immobile (e.g., the braiding operation is much faster than 
anyonic propagation), the contrast of the phase measurement 
only depends the probability that the loop l\ U h (or 1' 2 U 1' 4 ) 
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(see Fig. |4^) encloses an odd number of initial x-particles 
(or z-particles). An extra phase e' Q = — 1 will be accumu- 
lated from each loop satisfying this condition. Suppose the 
loop l\ U 1$ (or 1' 2 U 1' 4 ) encloses m faces (or m' vertices). If 
one pair of initial x-particles is uniformly distributed among 
jy2 x Qy2 _ 1 j 12 possible configurations, the probability for 

2m(N 2 — m) 

accumulating an extra phase is q m = n2 ^ N 2 [J f° r tne 1°°P 

l\ U I3. The probability reaches a maximum value sw 1/2 for 
m N 2 /2; meanwhile it vanishes for m = or A^ 2 , which is 
achieved by l\ — h- Therefore, the contrast for the fringes of 
(a<j)) v.s. will be reduced to 1 — p x (q m + q m i). 

If the anyons are highly diffusive (e.g., random anyonic 
propagation is very fast compared to the intervals between 
the control operations), we should avoid adding any anyons 
by applying string operators that act within the ground sub- 
space of the topological memory. As shown in Fig. [4};, we 
use generators of the encoded qubits associated with strings 
{Zzi,Zx2,/z3,/^ 4 } to implement the braiding operation. How- 
ever, any quenched anyons (if present) will quickly diffuse 
over the entire torus and completely wash away the fringe of 
Therefore, the remaining contrast with highly diffu- 
sive anyons is 1 — p. Note that imperfect string operators may 
also reduce the contrast, since they may introduce unwanted 
anyons to the topological memory with probability approxi- 
mately proportional to the length of the string. 

Extension to gauge theories 

This interferometric technique can be extended to measure 
abelian anyonic statistics for any Z t / gauge theory by intro- 
ducing the spin lattice Hamiltonian with d levels for each spin 
ll50l . One can still use a probe qubit to measure the statistical 
phase via controlled-string operations. The mutual statistical 
phase between a charge a £ Z f / and flux b E Z c i associated with 
the braiding operation is Z^X^Z^X^ = e a%ab / d . Here the 

string operator Z^ z (X b x ) is a product of Z" (X b ) operators of 
all the spins on the string Cz (Cx)> where Z and X are ele- 
ments of the generalized Pauli group. The operator Z can be 
implemented by phasing pairs of spin states at a time using the 
protocols in the main text for the appropriate duration at each 
stage. Equivalently, one can choose field polarizations and de- 
tunings so that only one of the d levels is strongly coupled to 
the cavity, then evolve that state for the appropriate time, and 
swap other states in, evolve, and swap out again. A total of 
d— 1 global gates suffice to simulate Z£ z . The same follows 
for the Xq operators, but one must perform a parallel Fourier 
transform operator F — YljeCxfj ^ rst on a ^ tne s P ms m me 
configuration, then implement Z b c then apply F~ 1 . 
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Supplementary Online Materials 

NOISE MODEL FOR TORIC-CODE HAMILTONIAN 

The toric-code Hamiltonian for spins on the edges ofNxN 
square lattice 



Hfopo — J y.-Aj J/.By> 



(8) 



is a sum of stabilizer operators A s — n/<=star(.s) 

crj and B p = 

n 7 -6boundary(p) <7 5 associated with the site (vertex) s and the 
plaquette (face) p, respectively. And the coupling strength 
J determines the energy gap between the ground and ex- 
cited states A ~ J, which is also the energy associated with 
the quasi-particle excitations. There are two types of quasi- 
particles: (1) z-particles that live on the vertices of the lat- 
tice and (2) x-particles that live on the plaquette. The quasi- 
particles do not change from one type to the other type, 
but there is a non-trivial topological phase associated with 
braiding of two quasi-particles of different type. We gener- 
ate and move these quasi-particles by applying string oper- 
ators, meanwhile during the interval between the string op- 
erators the quasi-particles will evolve under the toric-code 
Hamiltonian and various local perturbations from the envi- 
ronment. In the absence of perturbations from the environ- 
ment, the quasi-particles are immobile. However, the mo- 
bility of quasi-particles changes when we include local per- 
turbations, because the excited states with quasi-particles are 
highly degenerate and any small perturbation to the Hamil- 
tonian can change both the energy spectrum and the eigen- 
wavefunctions. 



Perturbation Hamiltonian 

In this section, we will consider a simple model that will 
induce diffusion of quasi-particles of the toric-code Hamil- 
tonian. We will consider the case that the perturbation is 
small compared to the energy gap A, so that he number of 
quasi-particles is still conserved. However, we will show that 
such small perturbation can lead to non-trivial dynamics in the 
manifold with fixed number of quasi-particles. The perturba- 
tion is described by the Hamiltonian 



Hpert — 



eGAll spins 



(9) 



where e is the label for spins that we sum over. The time 
dependent coefficient h e is the local field associated with spin 
e. 

We can also write the perturbation Hamiltonian by sum- 
ming over the plaquettes and their surrounding edges: 



Hpert — 2 E E hpftO^w 
P r|G5\£ 



where p is the label for the plaquette, r| is the label for the sur- 
rounding edges, 9i = {(0, 1) , (0, -1) , (1,0) , (-1,0)} is the 
set that includes four edges around the plaquette as illustrated 
in Fig. [6^. The combination of (p,T\) label the r|-edge of pla- 
quette p. We use hpft and op ^ to represent the local fluctu- 
ating field and the spin operator for edge (p,V[), respectively. 
Since each edge is shared by two neighboring plaquettes, both 
(p,T[) and (p + T[, — r|) represent the same edge. By definition, 

(s p+r\ ,-tv The pref actor 



we have h. 



l p+i\, 



-n and °£,n 



1 /2 in Eq. d9|) accounts for the double counting of edges. 



(a) 



-(o=i)- 



(-10) P (1,0) 

-to,-i)+- 




(10) 



FIG. 6: Relative position of edges with respect to (a) plaquette and 
(b) site. 



Effects from Perturbation 

We now consider the effects from the perturbation Hamil- 
tonian H pert , under the assumption that the number of quasi- 
particles is conserved. This assumption can be justified as 
long as the local fluctuating field has a small amplitude and 
low-frequency noise spectrum, compared with the energy gap 
from the topological Hamiltonian. For ground states of the 
toric-code Hamiltonian, there are no quasi-particle excita- 
tions. And the perturbation Hamiltonian can only induce vir- 
tual excitations of quasi-particles. Thus the leading non-trivial 
effect within the ground-state manifold is the Mh order per- 
turbative process, associated with hopping of virtual quasi- 
particles along a minimal non-contractible loop of length M 
Such Mh order process is suppressed by a factor {\h\/i±) N 
which decreases exponentially with respect to the system size. 
As expected, the ground states of the topological Hamiltonian 
should be robust against these local perturbations. 

The perturbation Hamiltonian acts very differently on the 
excited states. It can induce x-particles to hop to the neighbor- 
ing plaquettes, while having no influence to z-particles. We 
may define the creation operator of the x-particle at plaque- 
tte p as bp, which changes the stabilizer B p from +1 to — 1. 
Similarly, the annihilation operator b p changes B p from — 1 to 

+ 1. By definition, (£>],) 2 = {bp) 2 = 0. Applying o x to one 
of the edge of the plaquette will flip the sign of the stabilizer 
Bp, represented by b p + b^. Since each edge is shared by two 
neighboring plaquettes, cp ^ will flip both stabilizers B p and 
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Bp-pn which can be expressed as 



T m ~* ( b P + b l) ( b P+r\ 



(11) 



By conservation of quasi-particles, we can simplify the map- 
ping: 



b n b 



t 



blbn 



(12) 



Therefore, within the manifold of fixed quasi-particles, the 
perturbation Hamiltonian can be reduced to 



H'pert - Z hp,T\ {bp b P+\\ + b l b p+r\) j 



(13) 



where the local field h pX[ can also be interpreted as the tunnel- 
ing rate of quasi-particles from p to p + r|, or from p + T[ to 
P- 

We now use H ' n to study the dynamics of quasi-particles. 
Suppose at time t — we create an x-particle at plaquette p 



|<p ) = ^ |vac) 



And at a later time x the state becomes 



(14) 



l<Px) 

= exp[-iH pert x] b\\vac) 

1 - E ( h p^) 2 K \ vac ) + L (^Pfl 1 ) ^P+n l VflC ) 
n J i 



(15) 



where in the second equality we assume that h p r ^ is time- 
independent and expand only to the second order of h p ^i. 
The overlap between the initial and final states is 



(<po|<Px) « 1-£(A p ,t,t) 2 wexp 
1 



1 



(16) 



Therefore, for time-independent perturbation and relatively 
short waiting time (i.e., x <C h~ Jj), the probability for a quasi- 
particle to remain at the same position decreases quadratically 
with time. 

Generally, the local field h p ^ will depend on time and we 
should replace h PtTl X by J X h p ^ (t 1 ) dt' in the equations above. 
In particular, the local field h p ^ (f') can be a stochastic ran- 
dom variable. In the next section, we will study the dynamics 
associated with stochastic noise fields {h p ^\. 



Time-dependent Perturbation 

We now consider the local stochastic noise fields character- 
ized by the auto-correlation function 



where we assume that h pX[ is steady and the auto-correlation 
only depends on the time difference between the two sam- 
pling points. For simplicity, we will also assume independent 
local noise for different spins; that is the correlation function 
/ip.ri (0) h p i ]T |/ (t ) vanishes unless (p,r\) and (p',T|') represent 
the same spin. We can characterize the noise by using the 
power spectrum density, which is the Fourier transform of the 
auto-correlation function: 



For example, the Gaussian correlation function 

/(f)=^exp[- f 2 /x 2 ] 
has power spectrum density 



<=/, 



'■Mi c 



-exp[-co 2 /co 2 ] 



(18) 



(19) 



(20) 



where the typical amplitude of the stochastic field is the 
correlation time is x c , and the cut-off frequency is co c - = 2x~'. 
The power spectrum density vanishes for high frequency CO 3> 
co t .. 

We can estimate the probability for a quasi-particle to re- 
main in the same position, in the presence of local stochastic 
noise, 



pt = |(<po|<Pt)I 



h p ,r\ (t')dt' 



(21) 



We first evaluate the square average of the integral for the fluc- 
tuation field 



h pn (t')dt' 

fchpA ( f ')Vn (t")dt'dt h 
sin 2 (cox/2) - 



(cox/2) 
i2je/(o> = 0)t, 



-/(co)rfco 



(22) 



where the last step assumes that we are interested in a time 
scale much longer than the noise correlation time: x ^> X c . 
Plugging Eq. ( 22 1 into Eq. pi) , we get the probability 



p x w exp [-2zTx] , 



(23) 



where z= = 4 is the coordination number of the square 
lattice, and T = 2ti/(co = 0) is the diffusion rate to each 



neighboring position. For Gaussian correlation [Eq. (19l], 
the diffusion rate is 



f{t) = h PA (t')hp >n (t'+t), 



(17) 



r=2Vi^/co e 



(24) 



12 



Fringe Contrast for Interference Experiment 



We apply £/| operations at time x/2 and T. The final state will 
be 



We now consider how diffusion of quasi-particles affect the 
fringe contrast of the anyonic interferometry. The signal of the 
anyonic interferometry is attributed to the interference from 
path-ways with different braiding of anyonic quasi-particles, 
which is achieved by using controlled-string operations. Dur- 
ing the intervals between the controlled-string operations, the 
quasi-particles excited by the controlled-string operations will 
diffuse, and the final state of the topological memory will have 
components orthogonal to the initial state. Since the orthog- 
onal components do not contribute to the fringes, the contrast 
will be reduced. 

The anyonic interferometry is analogous to the Ramsey ex- 
periment in the following aspects. For both cases, we start 
with a superposition state |0) + 1 1) for some two level system. 
Then, we let the system evolve, and meanwhile it is also inter- 
acting with the environment (e.g., cavity mode and selected 
spins, or external magnetic field). Finally, we projectively 
measure the system in some basis |0) ±e ±,if 1 1). The reduction 
of the measurement signal is attributed to various decoherence 
processes. For Ramsey experiment, the dominant decoher- 
ence is induced by fluctuations of the external magnetic field 
and it is characterized by the dephasing time T 2 * . For anyonic 
interferometry, we can define a similar dephasing time 



2 zT 



(25) 



And the fringe contrast is equal to the averaged overlap func- 
tion, which can be expressed as 



|<cp |(pT>|«exp[-T/r 2 *]. 



(26) 



For stochastic noise with Gaussian correlation [Eq. ( 19 1], we 
have T* = 0) c / (2zV^l). 



Spin Echo Techniques 



qecho 



= ^ e -'-i7/ 2 ^(0^^ e -'/o V2 ^(0*'fot \ vac ) 



-il0 /2 Kr,('')-Hper,^/2+t')}dt' h f 



b p | vac) 



(29) 



For stochastic noise with Gaussian correlation [Eq. (|19|)], the 
averaged overlap function now becomes 



l(<Po|<P? 



echo\ 



'• exp 



(30) 



where T 2 ~ vV5a- Since T 2 ~ T 2 * x (t,^) 3/2 , for % c % h > 1 
we can extend the coherence time by using echo-techniques. 

Furthermore, we may introduce n pairs of time-reversal op- 
erations within the time interval T (e.g., the pulse sequence 



U; 



x/2n _ x/2n 



■Ui- 



). And the time dependence of the constrast 



Repeat n times 

becomes 



exp 



3 \ Tl 



The contrast reduction is further slow down by a factor of n 3 / 4 
and this is illustrated in Fig. 5 of the paper. 

General Perturbation Hamiltonian 

We now generalize the perturbation Hamiltonian by includ- 
ing stochastic fields along z direction 

H pert , g en= I K^e+K^e + K^e- (3D 
eGAll spins 

We introduce the creation (and annihilation) operators a] (and 
a s ) for z-particles at site s. Suppose the edge e connects two 
sites s and s + and it is also shared by two plaquettes p and 
p + T| (see Fig. 151. Within the manifold of fixed number of 
quasi-particles, we have the following mapping: 



Similar to NMR systems, we can also use spin-echo tech- 
niques to further suppress the stochastic noise. The essence 
of spin-echo is to apply an effective time-reversal operation 
in the mid of the evolution so that the noises from the two 
intervals cancel each other. For anyonic interferometry, the 
effective time-reversal operation for the perturbation H pert is 



n 

eGAll spins 



(27) 



This is because {H pert .U~} — and [H topo , U£\ = 0. For ex- 
ample, at time t = we create an x-particle at plaquette p 



m) 



■ b\ \vac) 



(28) 



- bpbp + n 



(32) 



(33) 



Thus, the o x e term leads to hopping of x-particles, while the 
a z e term leads to hopping of z-particles. The al term leads to 
hopping of dyonic particles (paired x-particle and z-particle 
sharing the same edge, e.g., a^bp) 



i Cm 



h h J 



(34) 

And the effective operator for <r e is still Hermitian. Note that 
a y e consists of terms like i (b p b 

p+r\ a s a l + ^ a s a s +i/ y p b P+r ^j, 
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and the minus sign for terms with different order is consistent 
with the phase associated with the braiding of anyons. 
The generalized effective Hamiltonian becomes 

perl, gen 



P T|G5\£ 



(35) 



E if % b P b l+t\+ b l b P+T\ > a s a ] 



eGAll spins 



where the first term induces hopping of x-particles and the 
second term for z-particles. These two terms commute with 
each other, except for the situation when (p,r[) and (s,Q rep- 
resents the same edge. The third term induces hopping of dy- 
onic particles. If we are studying diffusion property of quasi- 
particles that are far apart, there is essentially no dyonic parti- 
cles and we may safely neglect the effect from the third term. 
Therefore, as long as the diffusion does not induce braiding of 
quasi-particles, we can safely treat the diffusion for x-particles 
and z-particles as independent processes. 

We may also introduce the effective time-reversal opera- 
tion U£ = ILeAii spins for perturbations of £ e6A1] spins h z e G z e ., 
since they anti-commute {h z e o z e ,U£} — 0. We can combine U% 
and U% in a nested fashion to suppress the diffusion of both 
x-particles and z-particles: 



which can be further simplified as: 



ui- 



^ TT 



x/4 



x/4 



u x - 

" IT 



t/4 



(36) 



Note that since time reversal operations U% and also anti- 
commute with Y,e€ All spins h z e O z e , the nest combination of the 
two also suppress the diffusion of dyonic particles. Therefore, 
we are able to suppress diffusion induced by the general per- 



turbation Hamiltonian of Eq. ( 3 1 1 to higher order. 



Let us first consider the modification of TJ\, We refer to 
the boundary protruding edges in even rows as "even rough 
edges," and those in odd row as "odd rough edges." We define 

Un" = II eGAll spins ®l and U%° = II eeAU spins 

e^Odd rough edges e^Even rough edges 

so that U£' e and act on even and odd boundary edges, 
respectively. After this modification, both UY and E/£'° 
commute with all stabilizers, especially B p at the left and 
right "rough edges." Similarly, we label "even/odd smooth 
edges" associated the columns for the top and bottom "smooth 
edges," and modify U£ into U% e = U eeA a spins G z e and 

e^Odd smooth edges 

II eeAU spins cl, wmcn commute with all stabi- 

e^Even smooth edges 



U 1t 



lizers. 

We introduce the sub-sequence 



w 



for a, p = e or o. And finally the full echo sequence to sup- 
press anyonic diffusion for the surface-code Hamiltonian with 
boundaries is 

W e e (x/4) W e '° (x/4) W°' e (x/4) W°'° (x/4) . (37) 



Summary 

In summary, we have analyzed a simple noise model on top 
of the toric-code Hamiltonian. We have found that this noise 
model can explain diffusion of quasi-particles. For anyonic 
interferometry, the effect of quasi-particle diffusion is analo- 
gous to the dephasing of the Ramsey experiment. Based on 
this analogy, we have proposed a scheme to extend the spin- 
echo technique to the topological memory, which will further 
suppress the diffusion of quasi-particles. 



UNIVERSAL ROTATIONS ON THE TOPOLOGICAL 
MEMORY 



Time Reversal Operations for Surface-Code Hamiltonian with 
Boundaries 

We now consider the effective time-reversal operations for 
the surface-code Hamiltonian with boundaries. For the planar 
code on a square lattice (see Fig. la of the paper), at the left 
and right are "rough edges" where the stabilizer operator B p 
is a product of three a x spin operators associated with each 
boundary plaquette, while at the top and bottom are "smooth 
edges" where the stabilizer operator A s is a product of three 
a z spin operators associated with each boundary site. In con- 
trast to the stabilizers associated with interior sites and plaque- 
ttes, these boundary stabilizers anti-commute with the previ- 
ous echo unitary or U^. Thus, we have to modify the echo 
unitaries, so that they commute with all stabilizers. 



We can achieve universal rotations of the encoded qubit 
stored in the topological memory. For example, an arbitrary x- 
rotation Xq = e' 6x on the topological memory can be achieved 
via the gate teleportation circuit 









X 








j 



This circuit represents the following procedure: (1) use 
controlled-string operation A [X~\ to entangle the probe qubit 
(upper line) and the memory (lower line with a slash), (2) pro- 
jectively measure the probe qubit in a rotated basis, and (3) 
perform an encoded Pauli X gate over the topological mem- 
ory conditioned on the measurement outcome. Similarly, we 
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can also implement arbitrary z-rotation Zq = e lQZ on the topo- 
logical memory via the gate teleportation circuit 



|0> 



I*) -h 



^Z e 



H 



Since any rotation can be decomposed into a sequence of x- 
and z-rotations, the above two circuits suffice for universal 
rotations. 

The gate teleportation can be generalized to implement ar- 
bitrary unitaries generated by string operators. For string op- 
erator S, the unitary operation Sq = e' es can be achieved via 
the gate teleportation circuit 



1+) 



For the geometric phase gate scheme, we can actually im- 
plement rotations of the encoded qubit without the probe 
qubit, e.g., x-rotation of the encoded qubit can be decomposed 

as e M =D(-P)D(-ae'1 1 )D(P)D(ae ( 1 1 ) by choosing 

|aP|=8. 



